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A treatment of the characteristics of a plane wave propagating in a homogeneous, unbounded, 
and fully ionized plasma is given by employing the linearized, two-fluid continuum theory of plasma 
dynamics The plasma is assumed to be a macroscopieally neutral and loss-free mixture of gas of 
electrons and singly charged ions. The dispersion relations tor the general case oi propagation at 
an arbitrary angle to the direction of the static magnetic field are investigated without any limitation 
on the frequency and with emphasis on the coupling of the transverse and the longitudinal type waves. 
Simple analytical expressions for the dispersion relations in the various intervals of frequency and 
phase velocity, as well as for the frequencies where the transverse and the longitudinal type waves 
couple, are given in terms of the various parameters. This treatment does not cover the cases in which 
the propagation vector is in the close neighborhood of either the direction of the static magnetic field 
or that perpendicular to it. 



1. Introduction 

By far the largest number of applications to the study of the nature of wave propagation in 
the ionosphere are based on the magneto-ionic theory of Appleton and Hartree [Ratcliffe, 1959]. 
This theory is concerned with the treatment of the plane wave dispersion relations in a cold, homo- 
geneous and unbounded electron plasma. A number of extensions to the theory of Appleton and 
Hartree are available in the literature. The extension of the conventional magneto-ionic theory 
to the low frequencies by the inclusion of the motion of the heavy ions has been investigated by 
various authors [Astrom, 1950, 1951; Stix, 1957; Fejer, 1960; Booker, 1963; Story, 1956; Hines, 
1953, 1957, 1963; Seshadri, 1964a]. The magneto-ionic theory together with its hydromagnetic 
extension obtained by taking into account the motion of the heavy ions has been able to predict 
correctly the nature of the plane wave propagation in a plasma, under the circumstances in which 
the effects of the alternating electric field predominate over those of the pressure gradients. 

When the effects 6f the pressure gradients are not negligible, the magneto-ionic theory fails 
to provide an adequate description of the properties of wave motion in a plasma. Also when the 
pressures of the charged particles are taken into account, new types of wave motion are found to 
appear. There have been a number of investigations on the plane wave dispersion relations in 
a compressible plasma which is capable of sustaining not only the usual transverse electromagnetic 
type waves but also the longitudinal plasma waves. These investigations are all based on a linear- 
ized theory of plasma oscillations and the essentials of this theory have been summarized in a recent 
paper by Oster [I960]. The characteristics of a plane wave propagating in an unbounded single 
component warm plasma have been treated [Oster, 1960] only for the two special cases for which 
the propagation vector is either along or across the direction of the external magnetostatic field. 
Recently a systematic development of the theory of plane wave propagation in an unbounded, 
compressible electron plasma has been given for the general case in which the propagation vector 
makes an arbitrary angle with the direction of the magnetostatic field [Kieburtz, 1964; Seshadri, 
1964a and b]. 
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The plane wave dispersion relations in a compressible plasma wherein motions of both the 
electrons and the ions are taken into account have also been studied by Pai [1960, 1962], Watanabe 
[1961], and Tanenbaum [1961]. The treatments of these authors pertain to the case with no 
applied magnetic field or to simpler cases of propagation parallel or perpendicular to the direction 
of the external magnetic field. Various simpler aspects and special cases of the plane wave 
dispersion relations in a two component, warm plasma are contained in the books by Stix [1962], 
Allis, Buchsbaum, and Bers [1963] and Denisse and Delcroix [1963]. 

In this paper, the plane wave dispersion relations in a two component warm magnetoplasma 
are investigated in a systematic manner. The dispersion equation, which is essentially a quartic, 
is derived without making any approximations in the beginning and is cast in a form amenable 
to analysis by a perturbation procedure. As a natural consequence of the procedure employed 
in analyzing the dispersion relations, the coupling regions where the high phase velocity, pre- 
dominantly transverse waves merge with those of the low phase velocity, predominantly longi- 
tudinal waves become evident. Moreover, simple analytical expressions for various sections of 
the dispersion curves are obtained and this enables the proper understanding of the changes in 
the dispersion caused by the variations of the relevant parameters. 

It is only appropriate to mention the reasons for using the continuum theory in the present 
investigation. Basically, there are two methods of dealing with the plasma problems, namely 
the kinetic treatment based on the Boltzmann equation and the hydrodynamic treatment based 
on the various transport equations. The first method is more accurate but leads to greater mathe- 
matical complexities. The second method, though incapable of accounting for effects such as 
Landau damping, is more readily amenable to analysis. Moreover, the macroscopic treatment 
gives correctly the various waves and their phase velocities and this knowledge is bound to provide 
valuable guidelines for the more complete treatments based on the Boltzmann equation. 

2. Statement of the Problem 

Consider a homogeneous plasma of infinite extent. The plasma is idealized to be a lossless 
and macroscopically neutral mixture of gas of electrons and a single species of ions. The drift 
velocity of the electrons and the ions are assumed to be zero so that the plasma, as a whole, may 
be considered to be stationary. A uniform magnetic field P is assumed to be impressed exter- 
nally throughout the plasma in the z-direction, where x, y, and z form a right-handed rectangular 
coordinate system (fig. 1). It is proposed to restrict attention only to the linear, time-harmonic 
problem; the harmonic time dependence of the form e~ ia)t is implied for all the field components. 

Let M?, V e , and P e be respectively, the average number density, the velocity and the pressure 
of the electrons; and let TV,-, V\, and Pi be the corresponding quantities for the ions. Also let E 
and H be the alternating electric and magnetic fields. The linearized time-harmonic hydrody- 
namic equations of motion for the electrons and the ions are 

-ia)m e N e Ve = -N e eU?+l?eXzBo)-VPe (1) 

-i(omNiJ?i=Nie(E+vlx2Bo) - VP,- (2) 

where — e and m e are respectively the charge and the mass of an electron and -he and rm are the 
corresponding quantities for an ion. Since the plasma is assumed to be neutral, the number 
densities of the electrons and the ions are equal and hence, each may be set equal to TVo, for con- 
venience. The linearized equations of continuity combined with the equations of state are given by 

u 2 e m ( >N V • V e = io)P e (3) 

ufmiNoV •Vi = i(x)P i (4) 
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FIGURE 1. Geometry of the problem. 



where u v and m are the sound velocities in the electron gas and the ion gas respectively. In addi- 
tion, the electric and the magnetic fields satisfy the following time-harmonic Maxwell's equations 



V X E = ico/jLoH 



(5) 



S7xH=-ia>€oE+Noe(Vi-Ve) 



(6) 



where fio and eo are respectively the permeability and the dielectric constant of free space. 

It is desired to investigate the characteristics of a plane wave propagating in the plasma me- 
dium. Let k be the propagation vector such that zxk coincides with the y-axis and let be the 
angle that k makes with the direction of the static magnetic field, so that k. r = k sin 6 = kn, k y = 
and k z — k cos 6 = kl. Therefore, all the field components will have the spatial dependence of 
the form e ik{nx+lz) so that 



3 -i d -n d 
— =ikn: — — U; — : 

ox dy dz 



ikl. 



(7) 



On eliminating P e from (1) with the help of (3) and (7), an expression for V v can be obtained in terms 
of E. In a similar manner, the elimination of P\ from (2) with the help of (4) and (7) leads to an 
expression for V\ in terms of E. On substituting (5), (7), and the expressions for V e and V\ in terms 
of E into (6), the following set of equations specifying E x , E y , and E z is obtained. 



D n iD l2 

iDn D22 

Z>3i i D32 



013 




e' x 


1^23 




Ey 


Z)33J 




\E.j 



= [D][E]=0 



(8) 
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2 aiD; \ oj 2 ) 

jiCOci / k 2 ufl 2 \ 

CtiDA CO 2 ) 



k 2 lnCl co 2 np k 2 u 2 p ln o)L k 2 u 2 ln 

D 13 = D„ = -+ — + — — — — 

(o 2 co 2 a e D e (o 2 Q) 2 a\Di a) 2 

co 2 w 2 a e D e \ w 2 I w 2 aiDi\ w 

_ <OpeOice k 2 uUn (o'frwci k 2 uUn 



Ifc-^-l + ^(l-^ + ^l-2*3 (9) 



o 3 a e D e co 2 w 3 aiDi w 
k 2 n 2 C 2 oj 2 e ( k 2 u 2 n 2 \ co 2 ; / k 2 u 2 w 

wWe 

D n =1 ; — 1— I 2 — ; a„ = l =■ n = e,i. (10) 

Also ft>pf.(o)pi) and w ce (co r /) are the plasma and the gyromagnetic frequency of the electrons (ions) 
and are given by 

N e 2 eB 

(0 Dn = \ (Den — n — e,l 11) 



and Co = 1/ V/XoCo is the free space electromagnetic wave velocity. The determinant of [D] should 
be equal to zero in order that (8) may have a nontrivial solution and this condition leads to the 
dispersion equation. 

The two special cases of propagation along and across the static magnetic field are reviewed 
briefly before proceeding to the study of the dispersion relations for the arbitrary direction of 
propagation. 



3. Propagation Along the Direction of the Static Magnetic Field 

For the case of propagation along the direction of the static magnetic field, 1=1 and n^O. 
Therefore, €13 = €31 = €23 == €32 == -D13 == D31 = D23 = D32 = and €11=622, Du—Dzi and D33 = — €33. 
Consequently the general dispersion equation factors into the following three separate equations: 

Dn-Z)i2 = 0; D n + D 12 = 0; D 3 3=-€ 3 3 = 0. (12a, b, c) 

With the help of (9), (12a, b) may be shown to yield the following expressions: 

F=C (€, +€ 2 )- 1/2 : V=Co(e 1 -e 2 )- 112 (13a, b) 

wmere V — t is the phase velocity of the wave and 

_ 1 m 

61 " ~Ct 2 -R 2 ~n 2 ~R 2 m 2 (14a) 

R Rm 2 

62 n(n 2 -/? 2 ) si(n 2 -R 2 m 2 ) (14b) 



and 



(1 = — ' R = — ' m = — - (14c) 

(Ope (Ope m\ 
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In this investigation ra P /m, will be set equal to 1/1836 which corresponds to a hydrogen plasma. 
In general, m is very much small in comparison with unity. It is to be noted that fl is the nor- 
malized frequency and R is the normalized strength of the external magnetic held. The substi- 
tution of (14a) and (14b) in (13a) and (13b) yields 



EM(0) 



EM(e) 



= C 



-Co 



(Cl+R)(tl-Rm) 



ft 2 + (IR(1 - m) - (1 + m + R 2 m)_ 
(CL-R)(CL^Rm) 



n 2 -nR(l-m)-(l + rn + R 2 m)_ 



1/2 



1/2 



(15a) 



(15b) 



Note that in the phase velocities (15a) for the ordinary mode and (15b) for the extraordinary modes 
neither of the two acoustic velocities u e and u\ appear. Evidently therefore, (15a, b) pertain to 
the purely transverse electromagnetic modes with no coupling to the longitudinal plasma waves. 
These two modes are the same as in an incompressible plasma for the case of propagation along 
the static magnetic field and have been studied previously [Seshadri, 1964a]. In figures 2a and b 
the phase velocities given by (15a, b) are plotted for two values of /? 2 , namely R 2 = 1/2 and /? 2 = 4/3. 
In a similar manner, it may be easily shown with the help of (9) that (12c) yields 



(d 2 - 1 - m)V 4 - [fl 2 u 2 + mu 2 e (W - 2)] V 2 + Wmu\ = 0. 



(16) 



The dispersion equation (16) which does not contain the velocity Co of electromagnetic waves in 
free space obviously corresponds to the purely longitudinal plasma waves with no coupling to the 
purely transverse electromagnetic waves. Also (16) does not depend on the static magnetic field, 
as is to be expected since the purely longitudinal plasma waves are unaffected by the static mag- 
netic field for the case of propagation parallel to its direction. The study of the dispersion equation 
(16) has shown [Seshadri, 1965] that the electron plasma (EP) mode propagates for fl > 1 and its 
phase velocity approaches u ( in the limit of infinite frequency. On the other hand, the ion plasma 
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FIGURE 2a. Dispersion relations for propagation along FIGURE 2b. Dispersion relations for propagation along 

the static magnetic field (R 2 = h). , . . ~ . , / D> . 4\ 
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(IP) mode propagates for all frequencies and its phase velocity approaches u x in the limit of infinite 
frequency and the velocity V2 m, which is the velocity of the acoustic waves in the whole gas of 
charged particles, in the limit of zero frequency. In fact, the transition of the phase velocity of the 
ion plasma mode from the value V 2 u% to the value u% takes place rather rapidly near the frequency 
given by Q 2 = 2m. The phase velocities of the electron and the ion plasma modes, as given by 
the dispersion equation (16) are also plotted in figures 2a and 2b. The values of Co and u e are taken 
to be equal to 3 X 10 8 and 3 X 10 4 m/sec respectively. 



4. Propagation Across the Direction of the Static Magnetic Field 

For the case of propagation across the direction of the static magnetic field, / = and n = \. 
Therefore, as before, 613 = 631 = 623 — €32 = #13 = #31 =#23 — #32 = 0. Consequently the general 
dispersion equation factors into the following two separate equations: 

D 11 Z>22- J D? 2 = 0;Z>33 = 0. (17a, b) 

With the help of (9), (17b) may be shown to yield 

F EM (o) = C n(n 2 -l-m)- 1 / 2 . (18) 

Since neither of the two acoustic velocities, u e and u\ appear, (18) evidently corresponds to a purely 
transverse electromagnetic wave without any coupling to the longitudinal plasma waves. Also 
(18), which corresponds to the ordinary electromagnetic mode as indicated by the subscript EM(0) 
on V, is independent of the static magnetic field. As the ordinary electromagnetic mode has its 
electric vector in the direction of the static magnetic field, which exerts no force parallel to itself, 
its dispersion is unaffected by the static magnetic field for the case of propagation across its direc- 
tion. The phase velocity given by (18) is plotted in figure 3. 

With the help of (9) and (10), (17a) may be shown to yield the following' cubic equation in V 2 : 

A V 6 ^A 1 V 4 -hA 2 V 2 ^A 3 = (19) 

A = ft 4 - a 2 (2 + 2m + R 2 + R 2 m 2 ) + (1 + m + R 2 m) 2 (20a) 

A^-Cll^l^ 8 + m8)- W{1 + m + 8 + 4>m8 + m 2 8 + R 2 (l + m 2 

+ m8 + m 2 8)} + 2m8{l + m) + R 2 m(\ + m + 2m8 + R 2 m)] (20b) 

A 2 = £l 2 C 2 u 2 e [£l 2 (l + m + m8)-m(2 + 8 + m8 + R 2 + R 2 m)] (20c) 

A 3 = -WmuiCl (20d) 

when 

d = u 2 IC 2 . (20e) 

The dispersion equation (19) which contains Co, u e , and m, corresponds to modes that have trans- 
verse and longitudinal field components and in which the motion of both the electrons and the 
ions contribute to the wave propagation. It is possible to decompose the total field into three 
independent modes which have been designated as (i) the modified extraordinary electromagnetic 
mode [MEM(x)], (ii) the modified electron plasma mode [MEP], and (hi) the modified ion plasma 
mode [MIP]. The phase velocity of these three modes are given by (19), which is solved by a 
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FIGURE 3. Dispersion relations for propagation across 
the static magnetic field (R 2 = }). 



method similar to the one detailed in the next section and the following result is obtained: 



^MEM(J-) - Co 
MKEP = Co 



nw-i-/? 2 ) 



n 4 -nw+2)+ i 

ftW-i-/? 2 ) 



1/2 



for ft > n 2 



n 4 -n 2 (/? 2 H-2) + i 



1/2 



forn,<n<Vi+/? 2 



= aVCoZ 



1 1 1 . Q ^ -l-R 2 )} ^ 2 (ft 2 -!- /? 2 ) Co 



4fl 2 



2# 2 



W, 



1/2 



for III 2 — 1 — 7? 2 | 



Co 



= Ue£l(W ~ 1 - R 2 r m forVl+# 2 < ft < oo 

^Mip=Co[/? 2 m-ft 2 (l+/? 2 ] 1 / 2 forO<ft<ft« 

-^{ft 2 (l+/? 2 )-/? 2 m} + i [qift^l-f/? 2 )-/? 2 ^} 2 

- 4ft 2 C 2 w 2 {fl 2 - m(2 + /? 2 )}] 1/2 
m(2 + /? 2 )-ft 2 



= ftz^ 



[_ft 2 (i + /? 2 )-/? 2 m 
= {ft 2 (l+/? 2 )-/? 2 m}- 1 / 2 



12 



for ft a < ft < ft () 



ft 2 tf 2 



M^ 2 -™(2 + # 2 )} 



1/2 



(21) 
(22a) 
(22b) 



(22c) 
(23a) 



for|ft 2 -ft 2 |^tfi (23b) 
Co 



(23c) 



+ [ftX{^ 2 -^(2 + ^ 2 )} + 4mft 4 ^{ft 2 (l+/? 2 )-/? 2 m}] 1 / 2 

for|ft 2 -ft 2 |^m (23d) 

= ft V^^[ft 2 -m(2 + /? 2 )]- 1 / 2 for ft < ft < oo (2 3e) 
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where 



/ R 2 m \ ,/2 

ft « = TZ15I : fio=[/n(2 + /?2)]'/2. ( 24b, c) 



v l + 7? 
In obtaining (21) — (23), it has been assumed that 

ci»vi»m 



"0 



where V a — CoR vm is the Alfven wave velocity and U a — V 2m u e is the acoustic velocity in the 
entire gas of charged particles. The condition (25) is usually obtained im most practical cases 
and corresponds to the value of R 2 being of the order of unity. 

Although, as stated earlier, the three modes corresponding to (19) have both longitudinal 
and transverse field components, it may be easily verified that (i) the modified extraordinary elec- 
tromagnetic mode is predominantly transverse throughout its frequency range of propagation, 
2 < Q. < oo, (ii) the modified electron plas ma mo de is predominantly transverse for Oi< II < Vl4/? 2 
and is predominantly longitudinal for V 1 4- R 2 <fl< o°, and (hi) the modified ion plasma mode is 
predominantly transverse for < ft < ft a and is predominantly longitudinal for Cl a < ft < oo. 

In the limit of infinite frequency, the phase velocities of the MEM(jc), MEP and MIP modes 
approach Co, u e , and m respectively. Also in the limit of zero frequency, the phase velocity of 
the MIP mode approaches V a . The phase velocities of the three modes given by (21) — (23) are 
plotted in figure 3 for R 2 = 1/2. 



5. Propagation in an Arbitrary Direction With Respect to That of the Static 

Magnetic Field 

For arbitrary direction of propagation, the determinant of [D] does not factor into two simpler 
equations and this makes the analysis of the dispersion relations for the general case very com- 
plicated. With the help of (8), (9), and (10), the following quartic equation in V 2 is obtained: 

AoV 8 +A l V« + A 2 V 4 + A3V 2 +A 4 = (26) 

where 

A = (ft 2 ~ 1 - m) \i¥ - 2 (2 + 2m + R 2 + R 2 m 2 ) + (1 + 2m 4- 2R 2 m + m 2 + 2R 2 m 2 + R 4 m 2 )] (27 a ) 

A , = C 2 [ - H 6 (2 + 8 4- mS) 4- (1 4 {(4 + 4m + 26 + 6m8 + 2m 2 8) 

+ R 2 (2 + 2m 2 + 6/2 + m S + m 2 8 + m :i 8l 2 )} - (1 2 {(2 + 4m 

+ 2m 2 4- 8 4- lm8 + lm 2 8 + m 3 8) + R 2 ( 1 + 1 2 + 4m + 4m 2 

+ m?( 1 4- / 2 ) 4- 3m 67 2 4- m8 4- 4m 2 6 + m 3 6 4 3m 3 8/ 2 ) 

4 R 4 m 2 (2 4 8l 2 4- m8l 2 )} 4 2mS(l 4- 2m 4- m 2 ) ~h R 2 m( 1 4 I 2 

4 2m 4 2m/ 2 4 m 2 (l 4- 1 2 ) 4- 2m8 4- 2m6/ 2 4- 2m 2 8 

4 2m 2 8/ 2 ) 4 R 4 m 2 ( 1 4 1 2 4 m 4 ml 2 4- 2m8l 2 )]. (27b) 
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A 2 = Q [ft 6 (l + 28 + 2mb + mS 2 ) - ft 4 {(l + rra + 28 4- Sm8 + 2/n 2 8 

+ 2m8 2 -f 2m 2 8 2 ) + /? 2 (l+m 2 + 28/ 2 + 2m8 4 2m 2 8 + 2m 3 8P 

+ m8 2 / 2 + m*8H 2 )} + ft 2 {m8(4 + \m + 8 + 2m8 + m 2 8) 

+ K 2 (P + m + m 2 + m*P + 4m8/ 2 4- m8 + /ra8/ 2 + 4m 2 8 

+ m 3 8 + 5/n 3 8/ 2 4- 2w 2 8 2 / 2 + 2m 3 8 2 / 2 ) + R 4 m 2 (l + 28/ 2 

+ 2ro8P + m8 2 / 4 )} - # 2 / 2 m 2 {48(l + m) + J? 2 (l +/n+ 2m8(l + P))}]. (27c) 

A = - Cjaf [(1 6 (1 + m + 2m8) - fl 4 { 2m( 1 + 8 + /n8) + /? 2 (/ 2 + m + m 2 

+ Pin 3 + 2m8/ 2 + 2m 3 8/ 2 )} -f n 2 /{ 2 m/ 2 {2(l + m 2 + m8 

+ m 2 8) + R 2 m( 1 + m+ 2m8/ 2 )} - 2m 3 R 4 l\ (27d) 
^4 - Qu|n 2 m(fl 2 - # 2 / 2 ) (n 2 - R 2 /n 2 / 2 ) = ^ . (27c) 

The extremely lengthy and not altogether straightforward manipulations needed in obtaining (26) 
and (27) are omitted here for the sake of brevity. Note that 8= 10~ 8 . Since m, the ratio of the 
electron to the ion mass is of the order 10~ 4 , it is evident that 8 is of the order m 2 . With this 
knowledge it is clear that the coefficients of all the powers of d in (27) are in the form of a power 
series in m. It is seen from (26) and (27) that the dispersion equation (26) is a quartic in V 2 and a 
cubic in O 2 . Nevertheless it was found to be advantageous to solve (26) for V rather than for ft. 

It is proposed to solve (26) under the assumption (25). In almost all practical situations and 
particularly for the case of the ionosphere, (25) holds good and hence the following study is not 
restrictive from a practical point of view. In view of (25), it is legitimate to neglect m, 8, and R 2 m 
in comparison with unity and to omit 8 in comparison with R 2 . 

The dispersion as obtained from (26) is found to vary rapidly in the neighborhood of certain 
frequencies for propagation directions close to both along and across the static magnetic field. 
It is therefore desirable to study first the dispersion relations in the interior region which excludes 
both the axial (1 — I 2 > 8) and the transverse (/ 2 > m) boundary layers. The investigation of the 
dispersion relations for the axial and the transverse boundary layers will form the subject of a 
subsequent paper. 

The solution of the fairly complicated dispersion equation (26) is considerably simplified by 
the following perturbation procedure. If V of the order C or higher is only sought, it is found from 
(26) and (27), the leading terms are of the order C§ and the other terms are lower by a factor of at 
least 8. On retaining only the leading terms, the following simplified dispersion equation is 
obtained: 



where 



A«V*+AuV*+A»=Q (28) 

^oi=(ft 2 -i)[n 4 -n 2 (2+/? 2 )+i] 

- (ft 2 - 1) (O 2 - ft 2 ) (ft 2 - ft 2 ) (29a) 

/i 1? =-C 2 [2ft 6 -2ft 4 (2-f^ 2 )-hft 2 {2-h/? 2 (l + / 2 )}-« 2 m(H-/ 2 ] (29b) 
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An = cj [a« - n 4 (i + ft 2 ) + n 2 R 2 (i 2 + m )- rhw] 

= q(fi 2 - a 2 ,) (a 2 - ny (O 2 - ft 2 5i ) (29c) 






(30a) 



n|-/? 2 m 2 . (30b) 

It can be easily argued that (28) is valid only when |ft 2 — ftfj > 8 and |ft 2 — 0| f | > 8. Also, when 
ft is of the order m and V is of the order VraCo, the leading terms in (26) are of the order m 4 Cl and 
their retention gives (28) again. When (ft 2 — R 2 m 2 ) is of the order m 4 , the approximation (28) may 
be shown to be not valid. Consequently, (28) gives correctly the phase velocities of the order Vra 
Co when ft is of the order m or lower such that |ft 2 — R 2 m 2 \ > m 4 , and the phase velocities of the 
order Co for the higher values of ft such that |ft 2 — ft| f | > 8 and |ft 2 — ft|J > 8. 

The approximate dispersion equation (28) is the same as the one obtained from considerations 
of the magneto-ionic theory which neglects the electron and the ion pressures completely, and its 
analysis shows [Seshadri, 1964a] that there are two modes of propagation. The ordinary electro- 
magnetic mode propagates in the two frequency ranges ^ ft < ft 3 / and 1 ^ ft ^ oo and the extraor- 
dinary electromagnetic mode in the three frequency ranges, ^ ft < ft 5l , fti ^ ft < ft^ and 
ft 2 ^ ft ^ °°. The frequencies ft3 ? , Vi^u and fts* for which the phase velocity given by (28) goes 
to zero are called the resonant frequencies. From the foregoing discussion, it is clear that the 
magneto-ionic theory does not give correctly the dispersion near the resonant frequencies. The 
phase velocities specified by (28) become considerably smaller than Co in the close neighborhood 
of ft = ft 3 ; and ft = ft 4 j and very much smaller than VmCo in the close neighborhood of ft = ft 5 ;. 
In the close neighborhood of the resonant frequencies, the first term on the left side of (28) may be 
neglected in comparison with the other two terms with the result 

An (31) 

The solution of (28) gives 



^•'--zt^V*'-""**- (32) 

The upper and the lower signs in (32) correspond to the ordinary and the extraordinary electro- 
magnetic waves respectively. The phase velocities given by (32) and plotted in figure 4 for R 2 = 1/2 
are marked AB with various subscripts corresponding to the various branches. For the portions 
of the dispersion curves marked AB, the electron and the ion pressures play negligible roles as 
a result the corresponding field components are predominantly transverse in character similar to 
the pure electromagnetic wave. 

If V is of the order u e , the leading terms on the left side of (26) are of the order C%u\ and the 
remaining terms are lower by a factor of at least m. It follows, therefore, that for obtaining V of 
the order u e , only the leading terms need be retained resulting in the following expression for the 
phase velocity: 



V 2 — — — 

V ~ An (33) 
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fhere 



A :ii = -Ct ) uj[i¥ - ft 4 {2m + RHl 2 + m)} + 2mfl 2 R 2 l 2 - 2m 3 R 4 l 4 ] 



= (ft 2 - /? 2 / 2 ) (fl 2 - 2m) (H 2 - R*m 2 P). 



(34) 



It may be easily verified that (33) is valid only if |ft 2 — ft|J»8 and |ft 2 — ft| i |»8. Also it can be 
shown that (33) yields V of the order u e only for fl 2 >> m and |(1 2 — R 2 l 2 \>> m and these inequali- 
ties further limit the range of validity of (33). For ft 2 >> m, A 2 i, and A& may be simplified with 
the result (33) reduces to 



V 2 = 



u 2 mw-R 2 i 2 ) 



(n 2 



-Bid 



(" 2 -ftJi) 



(33a) 



which corresponds to a propagating mode only for ft 2 < il 2 < R 2 l 2 and fl 2 4 .< fl 2 < <*>. Note 
that flL 2 , </? 2 / 2 <D 4 2 . In figure 4, the sections of the dispersion curves specified by (33a) are 
marked CD with two different subscripts corresponding to the two different branches. Note that 
when ft tends to infinity, V s pecifi ed by (33a) approaches asymptotically the value a e . 

When Fis of the order wu e Co and Q 2 — fl 2 - or (I 2 — O 2 is of the order Vo, the leading terms 
in (26) are of the order ii*C% and the other terms are lower by a factor of Vo or smaller. The reten- 
tion of the leading terms alone yields: 



AuV 4 +A2iV*+A 3i = 0, 



(35) 



whose solutions are given by 



2A u ~ y i \2Au) Au' 



(36) 



It may be easily proved that Am/Ah < for ft in the close neighborhood of either (l& or 11 4/ . Hence 
only the upper sign in (36) will correspond to real values of V and consequently, the lower sign in 
(36) may be disregarded. 
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Figure 4. Dispersion relations for propagation at an 
arbitrary angle to the direction of the static magnetic 
field (R 2 = i, l 2 =£). 
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The phase velocities for ft 2 = ft| f and ft 2 = ft 2 ,, are obtained from (36) to be given by 



F ""f^- (37) 

W ith t he help of (29b) and (34), (37) is seen to yield real values for V of the order of magnitude 
vUeCo for ft 2 = ft 2 ,- and ft 2 = ft^ 2 - . Consequently, (37) with the upper sign, correctly gives the phase 
velocities when (ft 2 — ft 2 ? ) or (fj2_(|2.) j s f tne orc j er f magnitude Vo\ 

When (ft 2 -ft 2 ,) > V8~or (ft 2 -ft£) > Vs, it is found with the help of (29b, c) and (34) that 
Am/Au is of the order 8 lower than (J 2 /Mi/) 2 The square root in (36) may therefore be expanded. 
Also since for (ft 2 — ftf,-) ^> V8 or (ft 2 — ft 2 .) ^> V§, /4w and /4 2 i are of the same sign, (36) may be 

shown to yield 

When (ft 2 — ftf,-) <^ Vo or (ft 2 — ft 2 ,.) <^ Vo\ An and /<2? are of opposite signs and (36) may then be 
shown to reduce to 

^ = - A ~- (39) 

With the help of (39) and (31) the phase velocity specified by (36) is found to merge with that given 
by (32) for ft 2 — ft| ? . <^ Vo or ft 2 — ft| ? <^ Vo. In a similar manner, a comparison of (38) with 
(33) shows that the phase velocity given by (36) merges with that specified by (33) for ft 2 — ft 2 ,,- 
>Vd or ft 2 -ft 2 ,.^>\/o\ 

The phase velocities computed from (36) with the upper sign are also plotted in figure 4 
and the corresponding sections of the dispersion curves are marked BC with two different sub- 
scripts corresponding to the two different branches. Note that the curve marked A^BuCuDn forms 
one continuous dispersion curve. 

When V is of the order Vm u v , ft > m and |ft 2 — R 2 l 2 \ > m, the retention in (26) of only the 
leading terms, which are of the order CqUJtti, yields 



A 3 i 



V 2 =~ (40) 



When ft > m, (40) may be simplified with the help of (27e) and (34) to yield 

P = /m* 2 ft 2 /(ft 2 -2m). (41) 

Obviously (41) gives real values for V only for ft 2 > 2m. Further, when ft 2 — 2m is of the order m 2 , 
V given by (41) becomes of the order u € with the result the validity of (40) is ensured only under the 
additional restriction that ft 2 — 2m >m 2 . The phase velocities as computed from (40) are also 
plotted in figure 4. Note that (41) is independent of R. There are two branches corresponding 
to the following frequency ranges: (i) 2m < ft 2 < R 2 l 2 and (ii) R 2 l 2 < ft 2 < °° and these are respec- 
tively marked E$F$ and E4F4. 

Using arguments similar to above, it can be shown that for \fl 2 — R 2 l 2 | of the order of #i, V is 
of the order V m u ( > and is specified by 

AaV 4 +AxV*+A4i = 0. (42) 
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The solutions of (42) are given by 



SrV^J't- ^ 



In the close neighborhood of Q 2 = R 2 l 2 , it can be shown that 

A 3i _ u\ [il 4 - Vt 2 {2m + RH 2 + R 2 m} ± 2mR 2 l 2 ] 

A 2i Ct 2 (n 2 -R 2 ) (44a) 

A 4i _ muj(n 2 -R 2 l 2 ) 

A 2 ~ (W-R 2 ) ' (44b) 

When O 2 - # 2 / 2 is of the order m and R 2 > H 2 > # 2 / 2 , it is clear from (44b) that A 4i /A 2i < 0. Hence, 
only V+ given by (43) will be real. When Cl 2 — R 2 l 2 > m, A 4i /A 2i is smaller than {A^j2A 2i ) 2 by a 
factor of rn and the square root in (43) may be expanded. For fl 2 — R 2 l 2 > m and O 2 < R 2 , A 3i IA 2 i > 
and (43) may be then shown to reduce to 

i/2 = A 4 i_ mujft 2 

A:« (il 2 -2m)' (45) 

Comparison of (45) with (41) shows that V+ given by (43) merges with the phase velocity given by 
(41) for a 2 - RH 2 > m. When 2m < R 2 l 2 , Am/A* > and Ax/A* < 0. Hence V + and V- given by 
(43) are both real. If R 2 l 2 -W > m, A zi \A %i < and with the help of (43), it can be shown that 



Am 

~A 2i 



V2= - (46a) 



and 



_ A 4i = mu 2 n 2 

A 3i (il 2 -2m) ' ^ 46b ) 

An inspection of (33) and (46a) shows that F+ given by (43) merges with that given by (33) for 
R 2 l 2 > m. Similarly, it is seen that V- given by (43) merges with that given by (40) and (41) for 
R 2 l 2 -Ct 2 > m. Also, when tt 2 = R 2 l 2 , A 4i IA 2i = and A 3i /A 2i < 0. Hence, (43) yields 

V\=-~; V 2 _=0fortl 2 = R 2 l 2 . (47) 

A 2 i 

The phase velocities, V+ and V- given by (43) and plotted in figures 4 are marked D 4 E 4 and F5G5 
respectively. From the foregoing arguments, it is clear that the section of the dispersion curve 
marked D 4 E 4 merges continuously and smoothly with that marked C 4 D 4 on the low-frequency side 
and with the section £4/^4 on the high frequency side with the result the sections marked ^4^4, 
Z? 4 C4, C4D 4l D 4 E 4 and E 4 F 4 together form one smooth curve. Also the section F5G5 of the dispersion 
curve corresponding to V- given by (43) joins smoothly with the section marked E5F5 on the low- 
frequency side. Note that this curve goes to zero for £l = Rl. 

When 11 is of the order m, the phase velocity of the order vm u e may be shown to be specified 

by A 

^ — ^' (48) 



With the help of (29c) and (34), (48) may be simplified to yield for O of the order of m 

._ _ ,(n 2 -fi 2 / 2 m 2 ) 

V* = 2mu> e {n2 _ R2m2) - ( 49 ) 

591 



Only for < (1 < Rim and ft > Rm, does (49) yield real values for V. Note that V specified by 
(49) attains the value \2m u e l for fl tending to zero and goes to zero for (1 — Rim. It may be proved 
that (49) is valid even in the close neighborhood of £l = Rlm. When \Ct 2 — R 2 m 2 \ is of the order 
m\ V given by (49) becomes of the order u e and hence (49) is not valid. The phase velocity, as 
obtained from (49) is also plotted in figures 4 and the corresponding sections of the dispersion 
curves are marked by H\K\ and H5K5 respectively. 

When |fl 2 — R 2 m 2 \ is of the order m 3 , V is of the order u e and it may be shown to be specified by 

AuV 4 + A 2i V 2 + A : y, = (50) 

where 



A u = m(l + l 2 ) (51a) 



A 2 i = C 2 



2i — 1>5 



(51b) 



mi 2J rm)-R 2 l 2 m 2 
A,i = - 2C 2 u 2 l 2 m(tt 2 - R 2 l 2 m 2 ) . (51 c ) 

The two solutions of (50) are given by 



V 2 =-Al+ (AnY-43i 

Vl > 2 2A l r\\2A u ) Au (52) 

where V\ and V>> correspond to the upper and the lower signs in (52) respectively. For Cl 2 > R 2 mr, 
A-mIAu < and hence only \\ will be real. For ft 2 -R 2 m 2 > rrv\ (A 2i l2A u ) 2 > \Asi/Au\ and A\\\A%\ 
> 0; hence, (52) yields for V\ 

m A 3i n (£l 2 -RH 2 m 2 ) 

^=-A 2 r 2mU Uw-R*m*y W 

which is the same as that given by (48) and (49) showing that the phase velocity curve for V\ speci- 
fied by (52) merges with that given by (48) for il 2 — R 2 m 2 > m\ 

For RH 2 m 2 < ll 2 < R 2 m 2 , AsilAu < and hence again only V\ is real. Also if R 2 m 2 -tt 2 > m\ 
A\\\A-i\ < and therefore, (52) may be simplified to yield 

T/t) A21 

which is seen to be the same as (31). As a result, it follows that the phase velocity curve for V\ 
specified by (52) merges with that for V v given by (32). 

If il 2 < R 2 l 2 m 2 , AvIAu > 0, ItilAu < 0, and hence both V\ and V 2 are real. From (51) and 
(52), it is easily deduced that 

V 2 = -4 2 i 
' An 

and 

Clearly (55) shows that V>> specified by (52) is real and identical with that given by (48). The phase 
velocity curves computed from (52) are sketched in figure 4. The curve corresponding to V 2 is the 
same as the one marked H\K\ and that corresponding to V\ is marked BJ1$. From what has been 
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stated before, it is clear that the section B5H5 joins smoothly with the section A5B5 on the lower 
frequency side and with the section H5K5 on the higher frequency side. 

Following the same procedure employed before, it can be easily deduced that when O 2 is of 
the order m; V of the order \m u e is specified by the following equation: 

A2iV*+A 3i V*+A 4 i = 0. (56) 

For H 2 of the order m, the coefficients in (56) may be readily shown with the help of (29c), (34), and 
(27e) to be given by 

A~2i=l (57a) 

I, i = u 2 e (W-2m) ( 57b) 

A 4i = -muin 2 . (57c) 

The two solutions of (56) are 



n .=- 



3, 4 



Asi + // Am X 2 Ajj 
■JLT«-V\23CJ An (58) 



Since A\\\A^\ < 0, only V* corresponding to the upper sign in (58) is real. It follows from (58) that 

V\ = 2mu\ for il 2 <2m (59a) 



and 

mu 2 tt 2 



(it 2 -2m) 



for £l 2 >2m. (59b) 



It is seen from (41) and (59b) that the phase velocity curve for V :i given by (58) merges with that 
specified by (41) for (I 2 ^> 2m. For H 2 of the order m, (49) becomes identical to (59a) showing that 
the phase velocity curve for V3 specified by (58) for fl 2 < 2m merges with that given by (49). The 
phase velocity curve for F 3 , as obtained from (58) and plotted in figure 4 is marked K5E5. It is to 
be noted that the section of the dispersion curve marked K5E5 joins smoothly with the section 
H5K5 on the lower frequency side and with the section E5F5 on the higher frequency side. Also 
in the close neighborhood of £l 2=i 2m, V changes rather rapidly from the value V2ra u e which is 
equal to the acoustic velcoity in the gas of both the charged particles to the value vm u e which 
is equal to the acoustic velocity in the ion gas alone. This feature is the same as obtained for an 
isotropic two component plasma [Seshadri, 1965]. Further, it follows that the various sections 
of the dispersion curves marked respectively by A 5 B^ y BMh, H^K 5 , K 5 E$, E^F 5 , and F5G5 form one 
smooth curve. 

It is of interest to examine the dispersion for extremely low frequencies such that H 2 <^ R 2 m 2 . 
For this case, (26) may be simplified to yield 

P - R 2 mC 2 ( 1 + 1 2 ) V*. + R 4 l 2 m 2 C 4 V 4 ~ 2m*RH 4 C 4 u 2 V 2 - C 4 u 4 Ct 2 m 3 R 4 l 4 = 0. (60) 

Using arguments similar to those given previously, it follows from (60) that V of the order V. m Co 
is given by the following equation: 

V 4 - R 2 mC 2 (l + / 2 )P + RH 2 m 2 C 4 = 0. (61) 
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The solutions of (61) are 

Also the phase velocities of the order Vra u e are obtained from (60) to be given by 

V 2 = 2ml 2 u 2 e . (63) 

It follows therefore that in the low-frequency limit there are only three propagating modes with 
the phase velocities specified by (62) and (63). It is easily recognized that the first part of the 
section A4B4 of the dispersion curves in figure 4 corresponds to the fast hydromagnetic wave whose 
phase velocity is equal to the Alfven wave velocity V a = CoR\/m for all directions of propagation. 
The mode denoted by the first part of the section A5B5 of the dispersion curve in figure 4 is the slow 
hydromagnetic wave with a phase velocity which is equal to the Alfven wave velocity in the direc- 
tion of the static magnetic field #0 and which changes by a factor of cos 6 for the other directions 
of propagation. The mode denoted by H\K\ in figure 4 is the sound wave whose velocity in the 
direction of B is equal to the acoustic velocity V2/n u e in the entire gas of charged particles. 
The phase velocity of the sound wave, just like the slow hydromagnetic wave, changes by a factor 
of cos 6 for the other directions of propagation. The spatial dispersion of the three modes in the 
low frequency limit, as obtained here, is in accordance with the results obtained previously from 
magnetohydrodynamic considerations [Seshadri, 1964a]. 

The section of the dispersion curve denoted by KJZ$ for fl 2 considerably less than 2/n, has a 
phase velocity equal to the acoustic velocity in the entire gas of charged particles for all directions 
of propagation except in the transverse boundary layer. The static magnetic field has very little 
effect on the dispersion for this section of the phase velocity curve. 

Only the identification of the various modes given in figure 4 remain. In a two fluid compres- 
sible plasma there are only four independent modes of oscillation, and these are designated as 
follows: (i) the modified extraordinary electromagnetic mode MEM(Z), (ii) the modified ordinary 
electromagnetic mode, MEM(O), (iii) the modified electron plasma mode, MEP and (iv) the modified 
ion plasma mode, MIP. The branch A\B\ of the dispersion curve corresponds to the MEM(X) 
mode, since its cutoff depends on the strength B of the external magnetic field and since, in the 
limit of infinite frequency, its phase velocity approaches Co. The branch AJBi corresponds to the 
MEM(O) mode, since its cutoff is independent of B and since, in the limit of infinite frequency, the 
phase velocity of this mode also approaches Co- The branch A^B^C^D^ belongs to the MEP mode 
since, in the limit of infinite frequency, its phase velocity approaches the acoustic velocity u e in 
the electron gas. The branch A4B4C4D4E4F4 belongs to the MIP mode for, in the limit of infinite 
frequency, its phase velocity becomes equal to the acoustic velocity u\ in the ion gas. The remain- 
ing two branches of the dispersion curves do not extend to infinite frequency and are obviously 
parts of the above mentioned four modes. It may be easily argued that the branch A5G5 is a part 
of the MEM(X) mode. Also, it is appropriate to associate the branch H\K\ with the MEP mode. 

In an isotropic plasma, two out of the four possible modes are purely longitudinal and the other 
two are purely transverse. But in an anisotropic plasma, all the four modes contain both longi- 
tudinal and transverse field components, and for the purpose of distinguishing these from the cor- 
responding ones of the isotropic plasma, the word 'modified' is attached to their designation. 
Although the modes of anisotropic plasma, strictly speaking, contain both longitudinal and trans- 
verse field components, the procedure used in the analysis of the dispersion uncovers the fact that 
these modes are predominantly longitudinal in certain frequency ranges and predominantly trans- 
verse in the other ranges of frequency. The MEM(O) mode is predominantly transverse in its 
entire range of propagation, 1 < Q < 00. The MEM(Z) mode is predominantly transverse for 
< il < Rm and (I2 < H < °° and is predominantly longitudinal for Rm < Cl < RL The MEP mode 
is predominantly longitudinal for < fl < Rml, and Ct 4 i < H < 00 is predominantly transverse for 
Hi < fl < flu. The MIP mode is predominantly transverse in the frequency range < fl < il 3 i 
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and is predominantly longitudinal for CI-m < ft < o°. It is evident from the above discussion that 
in the neighborhood of the three frequencies ft = Rm, ft = ft.}, and = 124/, coupling of the trans- 
verse and the longitudinal type waves takes place. 



6. Concluding Remarks 

It is desirable to contrast the treatment of the plane wave propagation in a two component 
warm plasma contained in this paper with those aspects of the problem that are discussed in the 
books by Stix [1962], Allis, Buchsbaum, and Bers [1963] and Denisse and Delcroix [1963]. Apart 
from the careful consideration of the characteristics of plane wave propagation in a warm plasma 
through the use of Boltzmann equations, Stix [1962] has discussed the dispersion of the low fre- 
quency ion acoustic wave and the ion cyclotron wave. The low-frequency ion acoustic wave 
corresponds to the first half of the branch H\K\ of the dispersion curve, and the ion cyclotron wave 
consists of two parts, the first part corresponds to the branch Z? 5 // 5 and the second part to the 
second half of the branch H\K\ Since he treated only the low phase velocity and low-frequency 
approximation, Stix obtained a gap in the frequency spectrum in the range Rl < ft < R. Note 
however that in this frequency gap, the two Alfven waves propagate. 

A parametric representation was first introduced by Clemmow and Mullaly [1955] for ascer- 
taining the regions of propagation of the possible plane wave modes in an unbounded and homo- 
geneous magneto-ionic medium. The two parameters used by Clemmow and Mullaly [1955] are 
X = a) 2 ,/co 2 and Y 2 = a> 2 J(o 2 . In the analysis of the properties of wave propagation in a plasma, 
the determination of the propagation characteristics explicitly as a function of the wave frequency 
co is almost always desired. Such a determination is difficult if t he parameters X and Y 1 are used, 
since the wave frequency is mixed in both of them. Therefore, in the present analysis, instead 
of the conventional parameters, the following normalized frequencies are used: 

n = _a L;/?=5 > £ e. 

<*)p<> (Ope 

Allis, Buchsbaum, and Bers [1963] have considerably extended the Clemmow-Mullaly plots, for 
example, to apply to a two component plasma. The two parameters employed by Allis, Buchsbaum, 
and Bers are 

„ _ tope + w pi _ 9 0> C eft>ei 



CO" 



It is to be noted that a 2 is proportional to the number density and ft 2 is proportional to the square 
of the external static magnetic field. Note also that both the parameters a 2 and /8 2 , also contain 
the wave frequency. In order to overcome certain practical difficulties encountered in the con- 
struction of the graphical plots, Allis, Buchsbaum, and Bers [1963] use the assumed mass ratio 
m = m e /mi = l between an electron and an ion. Further, as the mass ratio vanishes, the param- 
eter /3 2 also vanishes. To avoid the above-mentioned practical difficulties, the graphical plots 
were constructed (but not included in this paper for the sake of brevity) in terms of the two para- 
meters ft 2 and R 2 A discussion of this new diagrammatic representation may be found in Seshadri 
[1964b] for a single component warm plasma and in Seshadri [1964a] for a two component cold 
plasma. This new parametric representation has two special features, namely, (i) the wave 
frequency co is contained only in one of the parameters and (ii) the mass ratio m does not appear 
in either of the parameters. The first feature facilitates greatly the study of the propagation char- 
acteristics explicitly as a function of frequency and the second feature enables the plots to be used 
with equal facility for both the single and the two component plasma, with the result it is possible to 
ascertain from the graphical plots the influence of the heavy ion motion on the propagation char- 
acteristics of a plane wave in a plasma. Moreover the plots in the ft 2 — R 2 parameter space can 
be employed as is for a multicomponent plasma. 
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Allis, Buchsbaum, and Bers [1963] have noted that the plane wave dispersion relation with 
electron and ion thermal motion included is sufficiently involved that no simple form has been found 
for it. Consequently, their discussion of the plane wave dispersion relations in a two component 
warm plasma proceeds by first neglecting the coupling of the electromagnetic and the plasma waves 
and the coupling being then described intuitively on the basis of the understanding of the corre- 
sponding coupling in a single component warm plasma. 

Denisse and Delcroix [1963] have considered various aspects of the characteristics of plane 
wave propagation in a two component warm plasma. For the case of propagation along and across 
the external magnetic field, they have plotted the dispersion curves in terms of x' and y' where 
x' = log [(co 2 , + w 2 f )/co 2 ] and y' = C 2 \V 2 for \C 2 /V 2 \ ^ 1 and y' = 1 + log (1 -flog |q/P|) for |C ( 2 } /P| ^ 1. 
In view of the different scales adopted for the different phase velocity regions and in view of the 
log log scale adopted for the low phase velocity region, a clear overall picture of the variation of 
the basic quantity, namely the phase velocity with the frequency, as is depicted in figures 2 and 3, 
is difficult to obtain from the work of Denisse and Delcroix. For the general case of propagation 
at an arbitrary angle to the direction of the static magnetic field, Denisse and Delcroix have at- 
tempted to extricate the rules governing the propagation of the coupled longitudinal and transverse 
modes. The dispersion relation has been expressed both as a quartic in u[ = (CfJV 2 — l) 2 ] and as a 
cubic in x[ = {a) 2 pe J c to 2 ,)/cD 2 ]. They have examined the vertical and the horizontal asymptotes cor- 
responding to x going to zero and infinity and thus have obtained respectively the limiting behavior 
of the dispersion for infinite and zero frequencies. They have not examined the dispersion over 
the whole frequency region nor are any dispersion diagrams for the arbitrary direction of propa- 
gation given. 

In contrast to the treatments contained in the books by Allis, Buchsbaum, and Bers [1963] 
and Denisse and Delcroix [1963], in this paper the dispersion equation for the arbitrary direction 
of propagation of the plane wave has been cast in a rather simplified form as given in (26) and (27). 
By systematically analyzing the dispersion equation in the various regimes of frequency and phase 
velocity, it has been possible not only to treat the coupling of the transverse and the longitudinal 
waves, but also to obtain simple analytical expressions for the dispersion relations in the various 
intervals of frequency and phase velocity, as well as for the frequencies where the transverse and 
the longitudinal type waves couple. Also, the dispersion curves are plotted in terms of the phase 
velocity as a function of frequency for a typical case in which the propagation vector makes an angle 
of 45 deg with the direction of the static magnetic field. 

In conclusion, it is appropriate to sum up that the plane wave dispersion relations in a two 
component warm plasma are worked out in this paper for the general case of propagation at an 
arbitrary angle to the direction of the static magnetic field. This treatment does not cover the 
transition regions corresponding to propagation vectors in the close neighborhood of either the 
direction of the magnetostatic field or that perpendicular to it. The treatment of the dispersion 
relations for the axial and the transverse boundary layers is reserved for a subsequent paper. 
Once the dispersion is known for all angles including the axial and the transverse transition regions, 
it will be possible to depict the spatial dispersion, that is, the variation of the phase velocity as 
a function of angle for various typical frequencies, as was done for the two component cold plasma 
[Seshadri, 1964a]. 
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